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Abstract 

We study the propagation of orientation waves in a director field with rotational in- 
ertia and potential energy given by the Oseen-Frank energy functional from the con- 
tinuum theory of nematic liquid crystals. There are two types of waves, which we call 
splay and twist waves. Weakly nonlinear splay waves are described by the quadrati- 
cally nonlinear Hunter-Saxton equation. Here, we show that weakly nonlinear twist 
waves are described by a new cubically nonlinear, completely integrable asymptotic 
equation. This equation provides a surprising representation of the Hunter-Saxton 
equation for u as an advection equation for v, where u xx = v x . There is an analogous 
representation of the Camassa-Holm equation. We use the asymptotic equation to 
analyze a one-dimensional initial value problem for the director-field equations with 
twist-wave initial data. 

Key words: Nonlinear hyperbolic waves, Liquid crystals, Variational principles, 
Integrable Hamiltonian PDEs. 
PACS: 04.30.Nk, 61.30.Cz. 



1 Introduction 



In this paper, we analyze a system of partial differential equations that models 
the propagation of orientation waves in a massive director field. The restoring 
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force for the waves is provided by the Oseen-Frank energy (2.1), used in the 
continuum theory of nematic liquid crystals, and the inertia is provided by 
the rotational inertia of the director field. 



The motion of the director field is then described by the variational principle 
(2.4). The corresponding Euler-Lagrange equation (2.5) is scale- invariant, and 
forms a non-dispersive hyperbolic system of wave equations. The system sup- 
ports two types of waves, which we call 'splay' and 'twist' waves, respectively. 
The principal nonlinear effects are that the wave speeds depend on the angle 
between the propagation direction and the director field, and that the twist 
waves generate splay waves. 

Splay waves were investigated by Saxton [1] and Hunter and Saxton [2]. In 
that case, the one-dimensional system of equations for the director field (see 
Section 4) reduces to a scalar wave equation for an angle <p(x,t), 



Here, the wave-speed a is a smooth, nonzero function of (p, given in (4.2), and 
the prime denotes a derivative with respect to ip. 

The wave equation (1.1) is obtained from the variational principle 



which is one of the simplest nonlinear generalizations of the variational princi- 
ple for the linear wave equation one can imagine. The effects of this nonlinear- 
ity include the formation of cusp-type singularities [3]. Smooth solutions can 
be extended by global weak solutions, but in sharp contrast with the more 
familiar case of entropy solutions of hyperbolic conservation laws [4], equa- 
tion (1.1) possesses conservative weak solutions that are compatible with its 
variational and Hamiltonian structure. The global well-posedness of the initial 
value problem for (1.1) for conservative weak solutions is proved in [5]. 

The use of weakly nonlinear asymptotics to study the behavior of solutions of 
(1.1) that consist of small, localized perturbations u(x,t) of a constant state 
fo leads (after an appropriate normalization, and provided that a'(<p ) ^ 0) 
to the Hunter-Saxton (HS) equation [2], 



ipu - a 2 (if)if x + a((f)a'(ip)<f 2 x = 0. 



(1.1) 




(iH + uu x ) x - -u 2 x = 0. 



(1.2) 
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This equation is derived from the variational principle 
t 1 

S J — \u x u t + uu x } dxdt = 0. (1.3) 



It is completely integrable [6,7,?], and possesses global dissipative and conser- 
vative weak solutions [9,10,11]. 

In this paper, we study the full system of director-field equations and show 
that qualitatively new phenomena arise which are not present for the scalar 
wave equation (1.1). The structure of the system is most easily seen in the 
case of one space-dimension; the director-field equations then reduce to two 
coupled wave equations (4.3)-(4.4) for angles ip(x,t), ip(x,t), given by the 
variational principle (4.1). The angle ip, with wave-speed a, corresponds to 
splay waves; the angle vp, with wave speed b, corresponds to twist waves. Both 
wave speeds a, b depend only on </?, and the wave equation for (p is forced by 
source terms that are proportional to quadratic functions of the derivatives 
of ip. If ip is constant, then the system (4.3)-(4.4) reduces to the scalar wave 
equation (1.1). 

Our main result is the following asymptotic PDE for a localized, weakly- 
nonlinear twist wave with amplitude v(x,t): 



(v t + uv x ) x = 0, (1.4) 

This cubically-nonlinear, scale-invariant system consists of an advection equa- 
tion for v in which the advection velocity u is reconstructed nonlocally and 
quadratically from v. It is obtained from the variational principle 

^ / - \v x v t + uvl + -u 2 x \ dxdt = 0. (1.5) 



The second dependent variable u(x,t) in (1.4) may be interpreted as the am- 
plitude of a splay wave that is generated nonlinearly by the twist wave, which 
then affects the velocity of the twist wave. 



Remarkably, the elimination of v from (1.4) implies that u satisfies the HS- 
equation 



(u t + uu x ) x - -u 2 x 



0. 



(1.6) 



Thus, a solution v of (1.4) is related to a solution u of (1.6) by v x = u xx . 
Under this change of variables, the variational principle (1.5) transforms into 
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a second variational principle for (1.6), distinct from (1.3) but compatible with 
it, that is associated with a Lie-Poisson structure (see Proposition 4). 



The HS-equation therefore arises from the system of director-field equations 
in two different asymptotic limits for two different types of waves, with two 
separate variational structures. The resulting bi-Hamiltonian structure ex- 
plains why (1.6) is completely integrable [6], and it follows that (1.4) is also 
completely integrable (see Section 6). Furthermore, as we show in Section 5, 
equation (1.4) may be solved explicitly by the method of characteristics; the 
key point is that the Jacobian of the transformation between spatial and char- 
acteristic coordinates satisfies an integrable Liouville equation. 

The correspondence between solutions of (1.4) and (1.6) is neither one-to-one 
or onto. Only convex solutions of (1.6) with u xx > can be obtained from 
solutions of (1.4), and solutions v of (1.4) for which v x has the same magnitude 
but different signs (which may depend on x) correspond to the same solution 
of (1.6). Moreover, because of the nonlinear nature of the transformation, dis- 
tributional solutions of (1.6) do not necessarily transform into distributional 
solutions of (1.4). Nevertheless, equation (1.4) provides an interesting repre- 
sentation of the HS-equation (1.6) for u as an advection equation for a new 
variable v. There is a similar representation of the closely related Camassa- 
Holm (CH) equation [12] (see (6.1)-(6.2) below, with Mu = u xx —u). 

The twist and splay waves exhibit a basic difference in their nonlinearity. Twist 
waves are linearly degenerate, meaning that the derivative of their wave speed 
with respect to the wave amplitude is identically zero, as stated in (2.12). 
This linear degeneracy may be seen in the asymptotic system (1.4), where 
the advection velocity u of v is independent of v, and in the one-dimensional 
system (4.3)-(4.4), where the wave speed b(ip) of ip is independent of ip. By 
contrast, excluding exceptional values of ip where a'(ip) = 0, splay waves are 
genuinely nonlinear, meaning that the derivative of their wave speed with 
respect to the wave amplitude is non-zero, as stated in (2.11). Despite the 
linear degeneracy of the twist waves, their interaction with splay waves leads 
to nontrivial, cubically nonlinear dynamics. 

Equation (1.4) differs from another, more obvious, cubically-nonlinear, scale- 
invariant modification of the HS-equation [2,13], 

(u t + u 2 u x ^j = uul, (1.7) 





given by the variational principle 
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Unlike (1.4), equation (1.7) is an asymptotic limit of the scalar wave equation 
(1.1), and arises when there is a loss of genuine nonlinearity at the unperturbed 
state ipo such that a'((p ) = but a"(ip ) ^ 0. In Section 3.3, we derive from 
the full director-field system a vector generalization of (1.7), given in (3.14) or 
(3. 16)— (3.17) , that describes the propagation of polarized orientation waves in 
the same direction as the unperturbed director field. 

We now outline the contents of this paper. In Section 2, we describe the 
system of PDEs for the director-field, and briefly compare it with some related 
variational field theories. We show that the system is hyperbolic, and study 
the linearized orientation waves. In Section 3, we summarize the asymptotic 
equations for weakly nonlinear splay and twist waves. 

In Section 4, we write out the one-dimensional director-field equations in 
terms of spherical polar angles, and construct an asymptotic solution of a 
one-dimensional initial value problem with initial data corresponding to a 
small-amplitude, compactly supported twist wave. This solution illustrates, in 
particular, the generation of splay waves by twist waves, and our main goal is 
to formulate equations for the resulting splay waves. 

There are two cases, depending on whether the twist waves are faster or slower 
than the splay waves (see Figure 2). When the twist waves are faster, they 
move at a constant velocity into an unperturbed director field, and we obtain a 
Cauchy problem for the splay- wave equation (1.1) with data for ip, ip x given on 
the space-like twist-wave trajectories. When the twist waves are slower, they 
are embedded inside the splay wave, and we obtain a free-boundary problem 
for (1.1). The twist- wave trajectories are not known a priori, and ip x satisfies 
a jump condition across the trajectories that follows from an integration of 
the asymptotic equation (1.4). 

In Section 5, we solve (1.4) by the method of characteristics, and use the result 
to solve the initial-boundary value problem for (1.4) that arises in Section 4. 
In Section 6, we show that (1.4) leads to the HS-equation, and describe its bi- 
Hamiltonian structure and integrability. We also consider the CH-equation. In 
Section 7, we derive an apparently non-integrable generalization of (1.4), given 
in (7.1)-(7.2), that applies to periodic twist waves with mean-field interactions. 



2 Director fields 

The system of wave equations we study here is motivated by the theory of ne- 
matic liquid crystals. We consider an unbounded anisotropic medium whose 
orientation (in three space dimensions) is described by a time-dependent di- 
rector field n(x, t) of unit vectors, n : M? x R — >■ § 2 . 
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For uniaxial nematic liquid crystals, the directions n and — n are identified, 
and then n : M 3 x R -> 1LP 2 . In most of the problems we consider here, the 
director field is a small perturbation of a given smooth field, and the global 
topological distinctions between S 2 and ILP 2 are not relevant. 

We suppose that the potential energy W of a director field n is given by the 
O seen- Frank energy functional [14] 

W (n, Vn) = \-a (div n) 2 + -f3 (n • curl n) 2 + ^7 |n x curl n| 2 . (2.1) 

Up to a null-Lagrangian, W is the most general quadratic function of Vn 
with coefficients depending on n that is invariant under the transformations 
n 1— > — n and 

x 1— > Rk, n 1— > i?n for all orthogonal maps i?. (2.2) 



The positive coefficients a, (3, 7 are elastic constants of splay, twist, and bend, 
respectively. For liquid crystals, we typically have < (3 < a < 7. 

In the one-constant approximation, a — f3 — 7, the Oseen-Frank energy re- 
duces (up to a null-Lagrangian) to the harmonic map energy, 

W(Vn) = i«|Vn| 2 . 



This energy does not depend explicitly on n, and is invariant under a larger 
group of transformations, namely 

x 1— > i?x, n 1— > for all orthogonal maps /?, 5. (2.3) 



The nonlinear effects we study here vanish in this case, and we assume that 
the elastic constant are distinct. 

We further suppose that the director field has rotational kinetic energy and 
is not subject to damping or dissipation. Its motion is then governed by the 
variational principle 

5 J j^n 2 - W(n, Vn)| d*dt = 0, n ■ n = 1, (2.4) 

R 3 xR 



where we have normalized the moment of inertia per unit volume of the di- 
rector field to one. 
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The Euler-Lagrange equation associated with (2.4) is 



n tt = aV(divn) — (3 {Acurln + curl (An)} 

+7{B x curln - curl (B x n)} + An, (2.5) 

where 

A = n • curl n, B = n x curl n. 



The Lagrange multiplier A(x, t) is chosen so that n ■ n = 1, which implies that 



A = — n t + a 



|Vn| — Icurlnl 



+ 2 



13 A 2 + 7B 2 + (a -7)divB. 



Equation (2.5) provides a natural, geometrical model for the propagation of 
orientation waves in an anisotropic medium, and is representative of a large 
class of variational systems of wave equations in which the wave speeds are 
functions of the dependent variables [15]. It is not applicable to standard liquid 
crystal hydrodynamics where the motion of the director field is dominated by 
viscosity and the effects of rotational inertia are negligible. Nevertheless, it is 
conceivable that the nonlinear phenomena analyzed here could be observed in 
high-frequency excitations of liquid crystals whose molecules possess a large 
moment of inertia (such as, perhaps, liquid crystals obtained from carbon 
nanotubes) . 

At a more general level, liquid crystalline phases are the result of a contin- 
uous breaking of rotational symmetry, and the orientation waves we analyze 
here may be regarded as associated Goldstone modes [16]. Similar phenomena 
should occur for non-dispersive Goldstone modes in other systems with con- 
tinuously broken symmetry in which the energy density associated with a set 
of order parameters \l/ a has the anisotropic form 

A{ ^ ) 



rather than the more commonly assumed isotropic Ginzburg-Landau form 
proportional to | X7\I> | 2 . For liquid crystals, the isotropic form arises only in the 
one-constant approximation. 

As liquid crystals illustrate, the Ginzburg-Landau energy density is not dic- 
tated by general symmetry arguments in anisotropic media. For long-wave 
variations, it is natural to retain the leading-order terms in the energy that 
are quadratic in the spatial derivatives of the order parameters, but the order 
parameters themselves may vary by a large amount, in which case one should 
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retain any dependence of the coefficients of the spatial derivatives on the order 
parameters. 

A similar situation occurs in classical field theories, where nonlinear sigma- 
models lead to wave equations with wave speeds that are independent of the 
dependent variables. This class of field theories includes the wave-map equa- 
tions from Minkowski space into § 2 to which (2.5) reduces in the one-constant 
approximation. On the other hand, general relativity leads to a form of non- 
linearity that is analogous to that of the general director-field equations: the 
wave operator in the Einstein equations acts on the metric and has coefficients 
that are functions of the metric. The effects of this nonlinearity in the Ein- 
stein equations are, however, much more degenerate than in the director-field 
equations [17]. 

We remark that the invariance properties of these equations also differ in an 
analogous way. Nonlinear sigma-models are invariant under separate transfor- 
mations of the independent and dependent variables, as in (2.3), whereas the 
gauge-invariance of the Einstein equations involves a simultaneous transfor- 
mation of the independent and dependent variables, as in (2.2). 



2.1 Orientation waves 



In this section, we show that (2.5) forms a hyperbolic system of PDEs, and 
describe the corresponding waves. We consider solutions of (2.5) of the form 

n(x, t) = n + n (x,t), 



where n' is a small perturbation of a constant director field no, linearize the 
resulting equations for n', and look for Fourier solutions of the linearized 
equations of the form 

n'(x,t) = iVe lk ' x -^n. 



Here, iV 6 R is an arbitrary amplitude, k G I 3 is the wavenumber vector, 
uj G K is the frequency, and n G M 3 is a normalized constant vector. 

We find that (see Appendix A.l) u satisfies the linearized dispersion relation 
u 2 - a 2 (k; no)] L 2 - b 2 {k; n )j = 0, (2.6) 



where, with k — |k 
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Fig. 1. Characteristic variety of (2.5). 



a 2 (k;n ) = a k 2 - (k • n ) 2 + 7 (k • n ) 2 , 
b 2 (k; n ) = (3 \k 2 - (k • n ) 2 l + 7 (k • n ) 2 . 



(2.7) 
(2.8) 



Thus, if a ^ /3, the characteristic variety of (2.5) consists of two nested el- 
liptical cones (see Figure 1). There is a loss of strict hyperbolicity when k is 
parallel to n , corresponding to a wave that propagates in the same direction 
as the unperturbed director field. 

The waves associated with the branch u 2 = a 2 (k;n ) carry perturbations of 
the director field in the same plane as n and k, in which n = R where 



We call these waves splay waves. We note that R is orthogonal to n , as 
required by the linearization of the constraint that n is a unit vector. 

The waves associated with the branch uj 2 = b 2 (k; n ) carry transverse pertur- 
bations of the director field orthogonal to n and k, in which n = S where 



We call these waves twist waves. 

There is a fundamental difference in the dependence of the splay and twist 
waves on n . From (2.7), we compute that the derivative of the splay-wave 



R (k; n ) = k - (k • n ) n . 



(2.9) 



S (k; n ) = k x n . 



(2.10) 
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speed in the direction of the perturbation carried by the wave is given by 



If a ^ 7, this quantity is nonzero provided that k is not parallel or orthogonal 
to n . On the other hand, from (2.8), we see that the derivative of the twist- 
wave speed in the direction of the wave is identically zero, 



By analogy with the terms introduced by Lax in the context of first-order 
hyperbolic systems of conservation laws [4], we say that a wave is genuinely 
nonlinear if the derivative of the wave speed in the direction of the wave is 
nonzero, and linearly degenerate if the derivative is identically zero. Thus, the 
splay waves are genuinely nonlinear when they do not propagate in directions 
parallel or orthogonal to the director field, and the twist waves are linearly 
degenerate. 



3 Weakly nonlinear waves 

In this section, we describe the asymptotic equations for weakly nonlinear 
splay and twist waves. We also consider the case of waves that propagate 
in the same direction as the unperturbed director field, when there is a loss 
of strict hyperbolicity and the orientation waves are polarized. The algebraic 
details of the derivations are summarized in Appendix A. 

3.1 Splay waves 

Weakly nonlinear asymptotics for genuinely nonlinear splay waves leads to the 
quadratically nonlinear Hunter-Saxton (HS) equation [2]. We summarize here 
the expansion for weakly nonlinear, non-planar splay waves. 

We look for a high-frequency asymptotic solution n e of (2.5) with phase 
$(x, t), depending on a small parameter e, of the form 




(2.11) 




(2.12) 




as e 



0+. 



(3.1) 
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In (3.2), n (x, t) is a smooth solution of (2.5) that is independent of e. 

We consider localized waves, such as pulses or fronts, rather than periodic 
waves. The resulting asymptotic solutions are valid near the wavefront $(x, t) = 
0, where 9 = 0(1), but they need not be uniformly valid in 9. Global solu- 
tions may be obtained by matching these 'inner' solutions for localized waves 
with suitable 'outer' solutions (see Section 4). For periodic waves, there are 
additional mean-field interactions (see Section 7). 

We define the local frequency oj(x,t) and wavenumber k(x, t) by 

w = -$ t , k = V$. (3.3) 



We assume that k ^ 0. It follows from the expansion that u, k satisfy the 
linearized dispersion relation (2.6), and we suppose that they satisfy the splay- 
wave dispersion relation, uj 2 = a 2 (k; no). The phase <£> then satisfies the lin- 
earized eikonal equation 

$ t 2 - a 2 (V$; n ) = 0. 



We assume that $(x, t) is single-valued and caustics do not arise. 



We find that (see Appendix A. 2) 
ni(0,x,t) = u(0,x,t) R(x,t), 



where R(x, t) is defined by (2.9). The scalar wave-amplitude function u(9, x, t) 
satisfies the HS-equation 

(u t + dL-Vu + Tuu e + Pu) e = )-Tu 2 e . (3.4) 



In this equation, a(x, t) is the linearized group velocity vector (a = Vk<^>), 
a = — [ak - (a - 7) (k • n ) n ] , 



r(x, t) is the genuine-nonlinearity coefficient (r = V no o; ■ R), 

a — 7 s 



r = - 



UJ 



(k-n ) k 2 -(k-n o y 
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and P(x, t) is given by 



P 



{ujR 2 ) t + div jujR 2 a) 
2uR 2 



where R 2 = k 2 — (k • n ) 2 . Equation (3.4) follows from the variational principle 



5 J -j^R 2 uq (ut + a.- Vu) + Tuu 2 e 



dOd^dt = 0. 



Introducing a derivative along the rays associated with $, 
d s = d t + Si-V, 

we may write (3.4) as an evolution equation for u along a ray, 

(u s + Yuu e + Pu) e = ]^u 2 e . (3.5) 

If J(x, t) is a non-zero ray-density function such that J t + div (Ja) = 0, then 
we have 

p i"R 2 )s Js 

2ujR 2 2 J' 
and the change of variables 



reduces (3.5) to (1.2). 



3.2 Twist waves 



We consider the propagation of weakly nonlinear twist waves with non-zero 
wavenumber vector k through an unperturbed director field n . We assume 
that a ^ f3 and k is not parallel to n . These conditions ensure the strict 
hyperbolicity of the system. 

As a result of their linear degeneracy, the effect of nonlinearity on the twist 
waves is cubic. We therefore look for an asymptotic solution of (2.5) of the 
form 
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n £ (x,*) = n [-±-L-L,t;e\, 



(3.6) 



n(0, x, i; e) = n (x, t) + e^n^O, x, t) + £ n 2 (0, x, t) + 0(e 3/2 ) (3.7) 

as e — > + , where n is a solution of (2.5). We assume that the local frequency 
and wavenumber (3.3) satisfy the linearized dispersion relation for twist waves, 
uj 2 = b 2 (k; n ), so that the phase $ satisfies the eikonal equation 

$ 2 - b 2 (V$; n ) = 0. 



We find that (see Appendix A. 3) 



ri! = vS, 

112 = - (^) ( k ' n o) uR + \ v 2 ( k x S ) > 



(3.8) 
(3.9) 



where R(x, £), S(x, t) are defined in (2.9), (2.10) and the scalar amplitude- 
functions u(9,x.,t), v(0,x.,t) satisfy 



(v t + b ■ Vv + Auvg + Qv) e = 0, 
u 9 e = v 2 . 

Here, b(x, t) is the group velocity vector, 



(3.10) 
(3.11) 



b=-[/?k-(/3- 7 ) (k-no)no], 

to 



A(x, t) is given by 

A= (^7)l 
cj(« — (3) 



k-n ) 2 A; 2 -(k-n o y 



and <5(x, t) is given by 

(cjS 2 ), + div (uS 2 b) 



Q 



2ujS 2 



where S 2 = k 2 — (k ■ n ) 2 . (With the normalization we adopt for R and S, we 
have R 2 = S 2 .) Equations (3.10)— (3.11) follow from the variational principle 



5 J ^loS 2 



v (v t + b- Vv) + A [ uv 2 + -u 2 e 



d6dxdt = 0. 
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The solution (3.6)-(3.9) consists of a leading-order twist wave with amplitude 
v and a higher-order forced splay wave with amplitude u that propagates 
at the twist-wave velocity. The amplitude and frequency of the forced splay 
wave are 0(e) and 0(l/e), respectively, which is the same scaling as in the 
weakly nonlinear solution for a free splay wave given in Section 3.1. The term 
proportional to v 2 in n 2 ensures that n is a unit vector up to the first order 
in e. 

The coefficient A of the nonlinear term in (3.10) is non-zero if (3 ^ 7 and k is 
not parallel or orthogonal to n . If k, n are constant and k is orthogonal to n , 
then there are exact large-amplitude traveling twist-wave solutions [18,19,20] 
and no weakly nonlinear effects arise. If k is parallel to n , then there is a 
loss of strict hyperbolicity and one obtains a system of asymptotic equations 
instead of a scalar equation. We consider this case in Section 3.3. 

Introducing a derivative along the rays associated with $, 
d s = d t + b • V, 

we may write (3.10)— (3.11) as 

(v s + Auv e + Qv) e = 0, u ee = v 2 . (3.12) 



If K(x, t) is a non-zero ray-density function such that K t + div (Kb) = 0, 
then we have 

V 2uoS 2 2K' 



and the change of variables 

u^——u, v^\——v, o e ^d x , d s ^——d t 
K V K ub 2 



reduces (3.12) to (1.4). 



3.3 Polarized waves 



The equations of motion (2.5) are invariant under spatial rotations and re- 
flections that leave n fixed. As a consequence of this invariance, there is a 
loss of strict hyperbolicity and genuine nonlinearity for waves that propagate 
in the same direction as n. The resulting polarized orientation waves are de- 
scribed by a cubically nonlinear, rotationally invariant asymptotic equation. 
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An analogous phenomenon occurs for rotationally invariant waves in first-order 
hyperbolic systems of conservation laws [21]. 



We suppose that the unperturbed director field n is constant, and consider 
waves with constant wavenumber vector k = kn parallel to n . We look for 
an asymptotic solution n e of (2.5) of the form 

, /k • x — uot \ 
n = n I ,x,t;e I , 

n(0, x, t; e) = n + e 1/2 ni(0, x, t) + en 2 (9, x, t) + 0(e 3/2 ) 
as e — > + . 

We find that J 2 = 7/c 2 (see Appendix A. 4), and 

n x = u where n • u = 0, 

n 2 = ~ (u • u) n . 

The leading-order perturbation u(9, x, t) satisfies the equation 

uj ,_, (a - 3)k 2 , . 
u et + -no • Vu e H — (u ■ u e ) d u 

+ (/9 ~J )A;2 {[(u • u) u e ] e - (u e ■ u e ) u} = 0, (3.13) 

which is derived from the variational principle 



■/Mu,.(u, + ^.Vu) + (^(u.„,f 



+ (g 7)^ ^ _ ^ ^ _ d0dxdt = 0. 
2cj I 



Making the change of variables 

d t + ^n • V i-> <9 t , <9 e ^ 9a,, 



and rescaling u, we can write (3.13) in a normalized form for u(x,t) G R 2 as 
Uxt + (/i - 1/) (u • u + 1/ [(u • u) u x ] -u(u x - u x ) u = 0, (3.14) 
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where 



a — 7 



li = 



v = 



7 



7 



The corresponding variational principle is 




Writing u = (u cosv , u sint>), we find that this variational principle becomes 



This result is consistent with what we obtain by expanding the one-dimensional 
variational principle (4.1)-(4.2) as ip — > 0, when 

a 2 ~ Oq + (a - 7)v? 2 , 6 2 ~ a 2 , + (/3 - 7)y? 2 , g 2 ~ y? 2 , 

and making a unidirectional approximation 9j ~ — ao^x in the resulting La- 
grangian. 

The Euler-Lagrange equations for (3.15) are 



This system is a coupled pair of wave equations for u and v. The radial mode 
u has velocity fiu 2 , so it is genuinely nonlinear when it / 0, while the an- 
gular mode v has velocity uu 2 , so it is linearly degenerate. If v is constant, 
corresponding to a plane-polarized wave, we recover the scalar cubic equation 
(1.7) for u. Nonlinear circularly polarized waves do not exist, however, since 
variations in the angular variable v force variations in the radial variable u. 



4 One-dimensional equations 

We consider a director field 

n(x,t) = (n 1 (x,t),n 2 (x,t),n 3 (x,t)) 




(3.15) 




u 



(3.16) 



(3.17) 
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that depends upon a single space variable x — x\. Writing 
n = (cos ip, sin p cos ip, sin ip sin ip) , 



where ip(x,t), ip(x,t) are spherical polar angles, we find that the variational 
principle (2.4) becomes 



*l\{<Pt- + K 2 " b\p)^] } dxdt = 0, 



(4.1) 



with 



a 2 ((f) = a sin 2 p + 7 cos 2 p, 

b 2 (p) = ( 5sin 2 < / 9 + 7cosV, (4.2) 
q 2 (p) = sin 2 <p. 



The Euler-Lagrange equation associated with (4.1) is a system of wave equa- 
tions, 



p tt - (a 2 p x ) x + aa'p 2 x + q 2 bb^ 2 x - qq' (V 2 - b 2 ^ 2 x ) =0, (4.3) 
At ~ (b 2 i; x ) x + ^ (p t ij t - b 2 p x ^j x ) = 0, (4.4) 

where the prime denotes a derivative with respect to p. The angle p corre- 
sponds to splay waves and the angle ip to twist waves. Both wave-speeds (a, b) 
are functions only of <p, and the wave equation for p> is forced by terms that 
are proportional to quadratic functions of derivatives of ip. 



4-1 The initial value problem 



Hunter and Saxton [2] construct an asymptotic solution p(x, t; e) of the initial 
value problem for the scalar wave equation (1.1) in —00 < x < 00 and t > 
with weakly nonlinear splay-wave initial data, 



p>(x, 0; e) = p + ef i^-j , p t (x, 0; e) = g i^-j , 

where /, g have compact support, and e is a small parameter. The solution 
consists of a superposition of right and left moving weakly nonlinear splay 
waves that originate from x — 0, whose width in x is of the order e. The splay 



17 



waves are separated by a slowly-varying, small-amplitude perturbation of the 
constant state ip that satisfies a linearized wave equation. 

Here, we construct an asymptotic solution for ip(x, t; e), ip(x, t; e) of the system 
(4.3)-(4.4) with weakly nonlinear twist-wave initial data, 

ip(x, 0; e) = (p , ip t (x, 0; e) = 0, 

1>(x, 0; e) = e^f , ^(x, 0; e) = e'^g , (4.5) 

where / and g are smooth, compactly supported functions. We suppose, with- 
out loss of generality, that the unperturbed constant value of ip is equal to 
zero. 

We use a 0-subscript on a function of ip to denote evaluation at ip = ipo. 
We assume that ao ^ bo, so that (4.3)-(4.4) is strictly hyperbolic at ipo. The 
system will then remain strictly hyperbolic, at least in some time-interval of 
the order one. We further assume that b' Q ^ 0, otherwise the leading-order 
nonlinear effects studied below vanish. In the case of the director-field wave 
speeds (4.2), these assumptions mean that ipo ^ nn/2 for jiGZ. 

Although ip is initially constant, it does not remain so. As we will show, the 
weakly nonlinear twist wave, whose initial energy 

oo oo 

- / {tf + b\M} dx = -J {g\6) + bim) d9 

— oo — oo 

is of the order one, generates a slowly-varying 'outer' splay-wave whose am- 
plitude is of the order one. 

We will construct an asymptotic solution of this initial value problem as fol- 
lows. 

(1) In a short initial layer, when t = 0(e), we use linearized theory. The 
initial data splits up into right and left moving twist waves. 

(2) For t = 0(1), nonlinear effects become important, and we use the method 
of matched asymptotic expansions, with different expansions for the twist 
and splay waves. 

(a) The twist waves are small-amplitude, localized waves, which vary 
on a spatial scale of the order e. We describe them by means of the 
weakly nonlinear asymptotic equations for twist waves derived above. 
We call these the 'inner' solutions. 

(b) Away from the twist waves, the leading-order solution is a large- 
amplitude splay wave, which varies on a spatial scale of the order 1. 
We call this the 'outer' solution. 
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Fig. 2. Characteristic structure for the solution of the initial value problem 
(4.3)-(4.5): (a) Fast twist waves; (b) Slow twist waves. The dashed lines are the 
trajectories of the twist waves, and the solid lines are the characteristic curves as- 
sociated with the splay wave. 

(c) We obtain jump conditions for the 'outer' splay-wave solution across 
the trajectories of the right and left moving twist waves by matching 
it with the 'inner' twist- wave solutions. 

(d) We obtain initial data for the nonlinear solution by matching it with 
the linearized solution as t — > + . 

The matching between the inner twist waves and the outer splay wave, and 
the structure of the nonlinear solution, depend on whether the twist waves 
are faster or slower than the splay waves (see Figure 2). In deriving the jump 
conditions for the splay wave below, we will consider these cases separately. 
For liquid crystals, we have f3 < a, meaning that twist deformations are not 
as 'stiff' as splay deformations, and the twist waves are slower. 

4-2 Initial layer 

First, we consider the solution of (4.3)-(4.4), (4.5) in a short initial layer when 
t = 0(e). Because of the finite propagation speed of the system, the solution 
is constant, with = <p , ip = 0, outside an interval of width of the order e 
containing x — 0. 

Near x = 0, we look for an asymptotic solution of the form 

V9 = $ (X,T) + O( £ ), ^ = s 1/2 $ 1 (X,T) + 0(e), X = -, T —-. 

e e 

Using this expansion in (4.3)-(4.5), we find that $ , ^i satisfy 
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%tt ~ (a 2 ($ )$ox) x + a($ )a'($o)$ox = 0, 

^itt ~ (b 2 (%)^ lx ) x + (%t^it - b 2 (%)% x q lx ) = 0, 

$o(X,0) =v^o, $ t(X,0) = 0, 

* 1 (X,0) = f(X), * 1T (X,0)=g(X). 

The initial value problem for $ has the constant solution $ = p , and 
therefore satisfies the linear wave equation 

^itt - tifoixx = 0, 

* 1 (X,0) = f(X), * 1T (X,0)=g(X). 

The solution is 

T) = F R (X - b T) +F L (X + b T) , (4.6) 
F R (9) = -f(9) + — Jg(Odt;, (4.7) 

1 1 °° 

F L {0) = -f{0)-— j g{Z)dH. (4.8) 
o ^ 

Since /, g have compact support, the functions F Rd , F L9 have compact support. 
4-3 Twist waves 

Next, we consider the propagation of the twist- waves for t = 0(1) through 
a possibly non-uniform splay- wave field ip(x,t). For definiteness, we consider 
the right-moving twist wave, which moves with velocity b(ip). The trajectory 
x — s n(t) of the wave satisfies 

ds R 

-dT = bR ' 

where an _R-subscript on a function of <p denotes evaluation at p = <PR(t), 
with (p R (t) = ip (s R (t),t). For the initial value problem, we have s R (0) = 
and p R (0) = (p . 

We introduce a stretched inner variable near this trajectory, 
= x - s R (t) 
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We find that the weakly nonlinear twist-wave solution of (4.3)-(4.4) is 



V = VR {t)+e V2 (9,t) + 0^' 2 ), 
where <p 2 , ipi satisfy 



(4.9) 
(4.10) 



2b R q R 



^269 



' QrPr^'i 
a R -b 2 ' ' 1 



'R 



tie- 



J R, 



o, 



(4.11) 
(4.12) 



We may transform (4.11)-(4.12) into (1.4) by a suitable change of variables, 
in which 9 corresponds to x, ipi to v, and <p 2 to u. 

Matching the twist- wave solution (4.10) as t — > + with the linearized solution 
(4.6) as T — > oo, we get the initial condition 



M0,0) = F R (9), 



(4.13) 



where F R is given in (4.7). Equations (4.11)-(4.12) are supplemented with 
suitable boundary conditions for ip 2 and ipi at 6 — ±oo, which we consider 
further below. 



The main result we need in order to obtain equations for the 'outer' splay 
wave solution is the following jump condition for ip 2 across the twist wave. 
We give a complete solution of the twist-wave equations in Section 5. 

Proposition 1 Suppose that ip 2 , ipi are smooth solutions of (4-ll)-(4-12) 
such that ipie(-,t) has compact support. Then 



d 
di 
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+ 



'a 2 - h 2 " 
a R °R 



R 



2 

¥26 



o, 



(4.14) 



where [■} denotes the jump from 9 = — oo to 9 = oo. 

Proof. Since ipig has compact support, it follows from (4.12) that ip 2 is a 
linear function of 9 for large negative and positive values of 9. Moreover, 



20 1 



R 



n 2 - h 2 
&R °Rj 



i,\ e d9. 



(4.15) 



21 



Multiplying (4.11) by ipw and using (4.12) to rewrite the result, we get 



m + 



bm 2g ra \ 2 

t- + = °- 



Integrating this equation with respect to 0, and using the fact that ipig has 
compact support, we find that 



d 7,o ,„ 1 f a% — b 2 K 
* + 2 IPT. 



I / < 



2 

¥>29 



b R Qr J 



Using (4.15) to eliminate the integral of ipf g from this equation, and rewriting 
the result, we get (4.14). □ 



For the coefficients in (4.2), the jump condition (4.14) becomes 

= 0. 



{b R tan (p R [ip 2e ] } + ^ (P - 7) sin 2 (p R 



¥29 



4-4 Matching: fast twist waves 



This case corresponds to < a < (3, when < a < b. Since the twist waves are 
faster than the splay waves, they propagate into a constant state (p = ipo, if) = 
ahead of them, and generate splay waves behind them. (See Figure 2(a).) It 
follows that the right and left moving twist waves move at a constant velocity 
along the trajectories x = b t and x = —b t, respectively. 

We consider the right-moving twist wave for definiteness. The appropriate 
inner variable is then 



e 



x — bot 



(4.16) 



The weakly nonlinear solution inside the twist wave must match for large 
positive 9 with the constant initial state ahead of the wave. This condition 
implies that (p 2 (9,t) and ipi(9,t) in (4.9)-(4.10), with ip R = ip , satisfy the 
boundary conditions 

V? 2 (oo, t) = 0, y?2e(oo, t) = 0, ^(00, t) = 0. (4.17) 

The derivative ip 2 e jumps from zero at 9 = 00 to a value 
<^2e(-oo,t) = a R (t) 
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at — — oo. It follows from the jump condition (4.14), with a R = ao, b R = bo, 
b' R = b' constants, that a R satisfies 

+ 2^ = °" 



From (4.12) and (4.13), we have ctr(O) = <tro where 

~i^ )_l F ™ de - 

The solution of this Riccati equation, 

a R (t) = am f - , (4.18) 



is defined for all t > 0, since ctro&o — when < a < b . 

In Section 5, we prove that when a < bo, equations (4.11)-(4.13), (4.17), with 
a R = a and so on, have a smooth solution defined for all t > 0. We note that 
the derivative ip 2 e(—oo,t) decays as t — > oo. This is a result of the fact that 
the twist wave radiates energy away from it in the form of splay waves. It also 
follows from the solution that the twist wave is a rarefaction, in the sense that 
its characteristics spread out with increasing time. 

A similar analysis applies to the left-moving twist wave, in which 
x + b t 

v = , 



and (p2o(0,t) = for 6 sufficiently large and negative. We find that <f 2 e = & L 
for 6 sufficiently large and positive, where 



with 



These 'inner' twist-wave solutions provide matching conditions for an 'outer' 
splay-wave solution ip(x,t). Using (4.9) with <p R = ipo to rewrite the condition 
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for the right-moving twist wave, 

V2o{6, t) ~ a R (t) as -> -oo 

where # is given by (4.16), in terms of the outer solution tp(x, t), and equating 
the outer limit of the inner solution with the inner limit of the outer solution, 
we find that 

ip x (x, t) ~ a R (t) asx^ b t~. 

Furthermore, the leading-order outer solution for (p is continuous across the 
twist-wave, and ip is higher-order in e. We obtain a condition for ip as x — > 
—bot + in an analogous way. 

Summarizing these results for the leading-order outer splay-wave solution 
tp(x,t), we find that ip = y?o is constant if x > b t or x < —b t. Inside 
the region — b t < x < b t, we find that <p satisfies (1.1), with data on the 
space-like lines x = ±b t given by 

ip(b t, t) = ip , <Px(b t, t) = a R (t), 

<p(-b t, t) = ipo, ip x (-b t, t) = a L (t), 

where a R , &l are given by (4.18), (4.19), respectively. 

Since the initial value problem for (1.1) is well-posed, this Cauchy problem is 
presumably solvable. The solution ip may form singularities, in which case it 
would have to be continued by a weak solution. 

4-5 Matching: slow twist waves 

This case corresponds to a > (3 > 0, when a > b > 0. Since the twist waves 
are slower than the splay waves, they generate splay waves both in front and 
behind them. As a result, the twist waves are embedded inside a splay- wave 
field. (See Figure 2(b).) The speeds of the twist waves depend on the splay- 
wave field, leading to a free-boundary problem for the trajectories of the twist 
waves, coupled with a wave equation for the splay wave that is subject to 
jump conditions across the twist-wave trajectories. 

We will not write out detailed asymptotic equations for the weakly nonlin- 
ear twist waves in this case, but we summarize the equations satisfied by the 
leading-order outer splay- wave solution ip(x,t). The main point is the deriva- 
tion of jump conditions for ip x across the twist-wave trajectories. 
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The right and left moving twist waves are located at x — sn(t) and x = si,(t), 
respectively, where 



^(t) = b(cp(s R (t),t)), ^jf(t) = -b(<p( 8L (t),t)), (4.20) 

s R (0) = 0, s L (0) = 0. (4.21) 

The solution is continuous across x = suit) and x — SL(t), so that 

Mil = °> Ml = 0. (4.22) 



Here, and below, we use [•]#, to denote the jumps across x = sn(i), 
x = s^it), respectively, meaning that 



Ma(*)= lir n - lim p(a;,t), 

[v]l(*)= Km p(a;,t)- lim v>(x,*)- 

a— >Si(t)+ x— >sz,(t) 



Considering the right-moving twist wave for definiteness, we have 



ip 26 (9, t) -> <r+(t) as -> oo, (4.23) 
¥29(0, t) -> <r_(f) as #^-00, (4.24) 

for some functions 0"+(i), <7-(i). The corresponding matching conditions for 

<f(x, t) are 

lim (fi x (x, t) — cr+(t), lim (p x (x, t) — cr-(t). (4-25) 

X-t-S+(t) s-»s-(t) 



From (4.14), (4.23)-(4.24), (4.25), and the analogous equations for the left- 
moving twist wave, we find that ip x satisfies the following jump conditions 
across the twist waves: 



-\ 


r f 4 - 




dt 1 






-\ 


r K - 




dt 1 





+ 



4 iJ 



ft 2 s 

°R 



R 



0, 



(4.26) 
(4.27) 



Furthermore, from (4.13) and (4.15), and their analogs for left-moving waves, 
we get the initial conditions 
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(4.28) 



Ml (0) 




(4.29) 



Finally, matching the outer solution with the initial, linearized solution, we 
find that (p(x,t) satisfies the initial conditions 



Summarizing, we find that the free-boundary problem for ip, sr, sl consists 
of (1.1) for (p(x,t) in — oo < x < oo, t > with the initial condition (4.30). 
The functions suit), siit) satisfy (4.20)-(4.21), and (p(x,t) satisfies the jump 
conditions (4.22), (4.26)-(4.29) across the curves x = s R (t), x = s L (t). 

We will not investigate this problem here. We remark, however, that Proposi- 
tion 3 in Section 5 implies that there is a smooth solution of the 'inner' asymp- 
totic equations for the weakly nonlinear twist wave (4.11)-(4.13), (4.23)-(4.24) 
whenever the derivatives cr±(t) = ip x (s^(t),t^j of the 'outer' splay- wave solu- 
tion of the free-boundary problem on either side of the twist wave are smooth 
functions of time. 



5 Method of characteristics 

In this section, we solve (1.4) by the method of characteristics. The explicit 
nature of this solution is related to the complete integrability of the equation, 
which is discussed in the next section. 

Proposition 2 Let (£, r) be characteristic coordinates for the PDE (1-4), 
where x = X(£,t), t = t, and write U(£,t) = u(X(£,t),t), V(£,t) = 
v (X(£, r), r). Then a formal solution of (1-4) is given by 



(p(x, 0) = (p , 



(p t (x,0) = 0, 



for i^O. 



(4.30) 



Tl 









(5.1) 



where A(£), B{r), F{£), G{r), H{r) are arbitrary functions. 
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Proof. Writing (1.4) in terms of characteristic coordinates (£, r) in which 
t — t and x T = u, we find that the PDE becomes 

X T = U, V 6r = 0, % - If-Uz = if, 



where J(£,t) is the Jacobian J = X^. It follows that V(£, t) = F(£) + G(r), 
where F, G are functions of integration, and 

J T = Ut, % - ^ = Ff. 



The elimination of U from these equations yields a PDE for J, 



Making the change of variables 77 = ??(£) where 77^ = F£, and J = —e K , we 
find that this PDE transforms into an integrable Liouville equation, 



KfjT c 



A" 



The general solution is 



where A(rj) and -B(r) are arbitrary functions. Integrating the equation X$ = 
—e~ K with respect to £ , we find that X is given by (5.1), which proves the 
result. □ 

Next, we consider (1.4) in —00 < x < 00, 



(v t + uv x ) x = 0, (5.2) 
u xx = vl, (5.3) 



supplemented with the initial condition 

v(x,0) = F(x), (5.4) 

and the boundary conditions 

u x (oo,t) = cr+(t), u x (-oo,t) = cr_(t). (5.5) 
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u,v 




Fig. 3. Schematic structure of the solution of the IB VP (5.2)— (5.5). 

We assume that F is a smooth function and that F x has compact support. 
Equation (5.2) then implies that v x (x,t) has compact support in x whenever 
a smooth solution exists, and (5.3) implies that u(x, t) is a linear function of x 
for sufficiently large positive and negative values of x. The boundary condition 
(5.5) specifies the corresponding values of u x . We illustrate the structure of 
the solution schematically in Figure 3. 

Computing the jump condition (4.14) for (5.2)-(5.3), we find that 



d 



— \u x ] + - 
dt [ J 2 



= 0, 



where [•] denotes the jump from x = — oo to x = oo. It follows that the data 
<t± must satisfy 

da + 1 o da- 1 9 „ N 
— r± + -al = — — + -a 2 . (5.6) 
dt 2 + dt 2 - v ; 



Moreover, integrating (5.3) with respect to x at £ = and using (5.4)-(5.5), 
we find that 



oo 



(7 + (0)- < t_(0)= J F 2 x dx. (5.7) 



The next proposition establishes the existence of smooth solutions of the IB VP 
(5.2)-(5.5) for compatible data F and a±. The solutions are not unique, since 
we may add an arbitrary function of time to v, and an arbitrary function of 
time to u (together with an appropriate time-dependent translation of the 
spatial coordinate x). We can remove this non-uniqueness by specifying, for 
example, u, v as functions of time at some value of x. 
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Proposition 3 Suppose that F : R — > R zs a smooth function and that F x has 
compact support. Also suppose that a+,a_ : [0, i*) — > R are smooth functions 
defined in some time interval < t < i*, where < £* < oo ; £/ia£ satisfy (5.6), 
(5.7). Then there is a smooth solution of (5.2)-(5.5), defined in — oo < rr < oo ; 

< t < u. 

Proof. If F(x) = F is constant, then <J+(t) = cr_(t), and a solution is u = 
xa+(t), v = Fq. We therefore assume that F is not constant. 

We then have from (5.7) that <r_(0) < cr + (0). Equation (5.6) implies that 



di 



(a+ - cr_) + - ((7+ + (7_) (cr+ - cr_) = 0, 



so <T-(t) < <r+(t) for all < t < t*. 



We choose constants ??_ < r] + < such that 



— oo 



and define the function r)(£) by 



— oo 
oo 



We also define a Jacobian J(£, r) by 



J = 



[(r? + -??)£+ + (77 -^)£_] 2 

(77+ - 77_) ((7+ - (7_) 



(5.8) 



where 




(5.10) 



(5.9) 
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One can verify that J = is obtained from (5.1) with 



B(r) = 



E + (t)-E-(t) 
r) + E + (T)-r)-E_(T) 



We then let 





u(£,t) = x t (s,t), v(0 = f(0- 

Since E-, E + > 0, < r) + , and rj- < r] < r)+, we see from (5.8) that J > 0. It 
follows that the transformation x = X(£, r), t — r between spatial and char- 
acteristic coordinates is smoothly invertible, and, according to Proposition 2, 
these expressions define a smooth solution of (5.2)-(5.3), as may be verified 
directly. We show that this solution satisfies the required initial and boundary 
conditions. 



First, at r = 0, we have E + = E_ = 1 and a + — o"_ = r] + 
(5.8) that J = 1 at r = 0, so x = £, and v(x, 0) = F(x). 



r]-. It follows from 



Second, using the equation 



_U L _Jr 

Ux ~ Xt ~ J ' 



we compute from (5.8) that 

^_ 2 (v+-v) d -t + (v-v-) d -t jK^z) , 

(V+ ~V)E+ + (V- V-) E- a+ - ct_ 



From (5.9)-(5.10), we have 



dE + 1 , , „ d£_ 1 



dr -i^--)^' ^T = 4 (c7+ - (7 - )E -' 



and from the jump condition (5.6), we have 
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Using these equations to eliminate r-derivatives from the expression for u x 
and simplifying the result, we get 

= (V+ ~ V) v~E+ + (V~ V-) v+E- 
(ri+-ri)E + + (ri-ri-)E- 



It follows that u x = er_ at x = — oo, when r] — r]-, and u x = a + at x = oo, 
when 77 = r^ + . □ 

For example, let us consider what happens when the derivative u x vanishes at 
x = — 00 or 1 = 00. If er_ = 0, then a + > satisfies the equation 

dcr+ 1 2 



which has a global smooth solution forward in time, 
o+(0) 



(7 + (*) 



l + (T+(0)*/2' 



It follows that (5.2)-(5.5) has a global smooth solution forward in time, which 
may be specified uniquely by the requirement that tt = u = 0ata; = -00. 
This case corresponds to the one that arises for the fast twist waves analyzed 
in Section 4.4. 

On the other hand, if a + = 0, then a_ < satisfies the equation 
da_ 1 



dt 2 
whose solution 

MO) 



l + <r_(0)f/2 



blows up as t t t*, where t* = — 2/cr_(0) > 0. Thus, a smooth solution of 
(5.2)-(5.5) exists only in the finite time-interval < t < t*. The derivative u x 
blows up simultaneously in the entire semi-infinite spatial interval to the left 
of the support of v x , so it does not appear possible to continue the smooth 
solution by any kind of distributional solution after the singularity forms. 
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6 Integrability and Hamiltonian structure 



In this section, we show that the twist-wave equation (1.4) is a completely 
integrable, bi-Hamiltonian PDE, and that if (u, v) satisfies (1.4) then u satisfies 
the HS-equation (1.6). 

We begin by describing the relation between (1.4) and the HS-equation (1.6). 
In order to describe the corresponding relation for the Camassa-Holm (CH) 
equation at the same time, we consider the following generalization of (1.4): 



(vt + uv x ) x = 0, (6.1) 
Mu = v 2 x , (6.2) 

where M is a self-adjoint linear operator acting on functions of x that com- 
mutes with d x . If M = d%, then (6.1)-(6.2) is (1.4). 

We suppose that u, v are smooth solutions of (6.1)-(6.2). Differentiating (6.2) 
with respect to t, using (6.1) to write v x t in terms of u, v and their spatial 
derivatives, then using (6.2) to eliminate v from the result, we find that u 
satisfies 

m t + mu x + (mu) x = with m = Mu. (6.3) 



If M = dl, then (6.3) is the HS-equation (1.6); if M 
the CH-equation [12], 



d 2 x - 1, then (6.3) is 



(ut + uu x ) x - l -u\ 



u t + 3uu 2 



Conversely, if m > is a smooth solution of (6.3) and v x = y/m, then v 
satisfies (6.1)-(6.2). 

Because of the nonlinearity of this transformation, difficulties may arise in its 
application to distributional solutions. For example, the function 

f if x < 

u(x,t) = < 

2x/t if x > 



is a weak solution of the HS-equation (1.6) in t > 0, and 
2 

u xx (x,t) = -5(x) 
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is non-negative in the sense of distributions. There is, however, no standard 
way to define a distribution v such that u xx = v x . 

The HS-equation (1.6) is bi-Hamiltonian and completely integrable [6], so 
(1.4) is also. Next, we consider the effect of the transformation v ^ u on the 
Hamiltonian structures of these equations. 

The system (6.1)-(6.2) is obtained from the variational principle 
5 J — j— v t v x — uvl + -wM-u j dxdt = 0. 



Variations with respect to u yield (6.2), and variations with respect to v yield 
(6.1). We may eliminate u by means of the constraint equation u = M _1 (v x ) 
to obtain a variational principle for v alone, 



6 J I {~ VtVx ~ 2^ M_1 (^)} dxdt = °- 



(6.4) 



The Euler-Lagrange equation for (6.4), 



Vt + M- 1 (vl)v x ] =0, 



(6.5) 



is equivalent to (6.1)-(6.2). Here, and below, we assume that operators such as 
M and d x are invertible; in the case of differential operators, this requires the 
addition of suitable boundary conditions which we do not specify explicitly. 

Making a Legendre transform of the Lagrangian in (6.4), we get the corre- 
sponding Hamiltonian form of (6.5), 

vt = 9- 1 {jfj , H = \j vlM^ (vl) dx, (6.6) 



where d x 1 is the constant Hamiltonian operator associated canonically with 
the variational principle (6.4). 

Proposition 4 Let denote the Poisson bracket of junctionals T , Q of 

v associated with the constant Hamiltonian operator d~ l , 



Under the change of variables u — M 1 iyl), where M is a self-adjoint linear 
operator, the bracket (6.7) transforms formally into a Lie-Poisson bracket 
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J = -2M- 1 (md x + d x m) M~\ 
where m = Mu. 



(6.8) 
(6.9) 



Proof. First, we consider the nonlinear change of variables m — v 2 x . Variations 
of v of the form v £ = v + ek + . . . lead to variations m e = m + eh + . . . of m 
where h = 2v x k x . For any functional T of m, with T £ = jF(m £ ), we have 



d_ 

de 



£ = 



/ 



5m 



hdx. 



Writing h in terms of k and using the skew-adjointness of d x , we compute that 



d_ 

de 



£=0 \ fX 



Since 



d_ 

de 



11^ 



£ = 



we conclude that 

— = -2 (— vm 

Sv \ 5m 



Using this equation in (6.7) and integrating by parts, we get 



) 5m 



= — 2 / (md x + d x m) —— dx. 
J 5m 5m 



Making the linear change of variables u = M l m in this expression, and using 
the self-adjointness of M, we get (6.8)-(6.9). □ 

The Hamiltonian form of (6.3) for u corresponding to (6.6) for v is therefore 
'&H\ 



u t = J 



du 



TC — - [ uMu dx, 
4 J 
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where J is given in (6.9). 

To give a second Hamiltonian structure for (6.1)-(6.2), we define a skew- 
adjoint operator K, depending on v, by 

K = v x d~ x M~ x v x , (6.10) 

where M is a self-adjoint linear operator commuting with d x , as before. 
We find that the operator (6.10) satisfies the Jacobi identity if the quantity 
(gh x - hg x ) Mf x + (hf x - fh x ) Mg x + (fg x - gf x ) Mh x 

is an exact x-derivative for arbitrary functions /, g, h. This condition holds 
for M — d 2 , since 

(dh x hgx) f xxx \{jjh x hg x ) f xx ] x ~\~ hf xx g xx gh xx f xx , 

and the terms that are not exact derivatives cancel under a cyclic summa- 
tion. The condition also holds for M = d 2 — 1. Moreover, in those cases, 
{c\d~ l + C2K) satisfies the Jacobi identity for arbitrary real constants ci, C2, 
so that d' 1 and K define compatible Hamiltonian structures. 

The Hamiltonian form of (6.5) with respect to K is 

v t = K(j?J, V = Jv 2 x dx. (6.11) 

If M = d 2 xl then K = v x d~ 3 v x , and (6.11) is equivalent to (1.4). 
Under the transformation m = v 2 ., equation (6.11) becomes 




K = — (md x + d x m) (d x 1 M (md x + d x m) , 
which gives (6.3). 

When M = d x , equations (6.6) and (6.11) provide a bi-Hamiltonian structure 
for (1.4). One can then obtain an infinite sequence of commuting Hamiltonian 
flows by recursion. We will not write them out explicitly here, but we remark 
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that among them is a Hamiltonian structure for v which maps to the Hamil- 
tonian structure for u canonically associated with the variational principle in 
(1.3). 



Lax pairs for (1.4) follow directly by transformation of the Lax pairs for the 
HS-equation [6]. We define 



Then, using the identity fd x 1 — d x 1 f = d x 1 f x d x 1 , we compute that the Lax 
equation L t = [L, A] is equivalent to 



which may be rewritten as (1.4). Alternatively, we can set u = d x 2 (v 2 .) in the 
original Lax pair. 



7 Periodic twist waves 



Spatially periodic splay waves are described by the following version of the 
HS-equation [22,23] 



Here, u(x, t) is a periodic function of x, and angular brackets denote an average 
over a period. The wave also drives a mean-field, which evolves on the same 
time-scale, t — 0(1), as the wave. 

The interaction between a weakly nonlinear twist wave and a mean-field is 
more complicated because the mean-field evolves on a faster time-scale than 
the wave. This is a consequence of the fact that the nonlinear self-interaction 
of a weakly nonlinear twist wave is cubic, but the mean-field is driven by 
quadratic nonlinearities. 

In this section, we derive the following generalization of the twist- wave asymp- 
totic equation (1.4) that applies to periodic waves: 








(v t + uv x ) x + n{v 2 x )v = 0, 
u xx = v 2 x - (v 2 x ). 
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Here, u(x, t), v(x, t) are periodic functions of x, which we assume to have zero 
mean without loss generality, and /i is a constant that cannot be removed by 
rescaling. We will derive (7.1)-(7.2) from the one-dimensional wave equations 
(4.3)-(4.4), but a similar derivation would apply to more general systems. 

It is interesting to note that the mean-field interaction introduces a dispersive 
term of Klein-Gordon type into the evolution equation (7.1) for v, despite 
the fact that the original system is scale-invariant and non-dispersive. The 
scale-invariance is preserved by the fact that the coefficient of the dispersive 
term is proportional to the mean wave-energy (or momentum) (v%), which is 
a constant conserved quantity. 

The mean-terms prevent the elimination of v from the system (7.1)-(7.2) by 
cross-differentiation, as is possible in the case of (6.1)-(6.2). Instead, one finds 
that 



(iH + uu x ) x - -ul - (v 2 x ) (2u + fiv 2 



0, 



which suggests that (7.1)-(7.2) is not completely integrable when /i ^ 0. 



7. 1 Derivation of the periodic equation 



We consider spatially-periodic solutions of (4.3)-(4.4) with period of the order 
e in which ip has amplitude of the order e 1 / 2 , where e is a small parameter. 
The corresponding time-scale for the nonlinear evolution of the xjj-wsve is of 
the order 1. As we will see, the -^-wave generates a mean <y?-field which evolves 
on a time-scale of the order e 1 ! 2 . This mean field modulates the speed of the 
ip-wme on the same time-scale. 

We therefore introduce multiple-scale variables 



x - e l / 2 s (t/e 1 / 2 ) 
e 

t = t, 

where s(r) is a suitable phase function, and look for an asymptotic solution 
of (4.3)-(4.4) of the form 

ip = ip (r)+e^2 (6,r,t) + 0(e 3/2 ), 

i> = e 1 ' 2 ^ (6, r, t) + e^ 2 (6, r, t) + 0(e 3 / 2 ), 
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where all the terms are periodic functions of 9. We will also require below that 
the terms are periodic functions of r. 

We use this ansatz in (4.3)-(4.4), expand derivatives as 
1 1 1 

9 t -> S T dg + — yrdr + 8 t , 8 X -> -8g, 



Taylor expand the result with respect to e, and equate coefficients of powers 
of e 1 / 2 to zero. 

We consider first the ^-equation (4.4). At the order e~ 3//2 , we obtain that 
(4 - b o) iPiee = 0, 

where the 0-subscript on a function of (p denotes evaluation at tp — ipo. It 
follows from this equation that s 2 = b$. For definiteness, we consider a right- 
moving wave and assume that 

Sr = b , (7.4) 

where b > 0. 

At the order e -1 , we obtain that 

2s T V>i0r + ( s TT + -^s T ] ijj ie = 0. (7.5) 

We may choose ip\ so that it is a zero-mean periodic function of 9. It follows 
from (7.4) and (7.5) that 



Qo(r)bl /2 (r) 



M°,r,t)= Zll , ( 7 - 6 ) 



where v(9,t) is a zero-mean periodic function of 9 which is independent of r. 
At the order e~ 1//2 , we obtain that 

2s r ?/W + (s TT + -^St^ ^28 + 2s T ipi0t + (2b b' ip 2 ilj ie ) g - ip lTT = 0. 

Using (7.4), we may write this equation as 

(qobl /2 ip 2 e) T + qobl /2 ipwt + (Wo^Or^ifl) fl - -^r^^irr = °- ( 7 - 7 ) 



38 



We will return to (7.7) after we expand the (^-equation. 

The leading-order terms in the expansion of the (^-equation (4.3) are of the 
order e^ 1 , and give 

(s 2 r - do) ¥268 + qlhb'Q^le + Vott = 0. 



Using (7.4) and (7.6) in this equation, we get 
(b 2 - afj ip 2 ee + b' v 2 + y?orr = 0. 



(7.8) 



Averaging this equation with respect to 9, we find that 
V^orr + (ve)b'o = 0, 



(7.9) 



where the angular brackets denote an average with respect to 0. As we will 
see, for smooth solutions, the quantity (vj) is a constant independent of t, so 
(7.9) is consistent with the ansatz that ipo depends only on r. 

Equation (7.9) provides an ODE for (p , corresponding to motion in a potential 
proportional to the twist- wave speed b = b((p ). For definiteness, we assume 
that the solution of (7.9) for <po(r) is a periodic function of r. We then require 
that all other terms in the expansion are periodic functions of r. 

Subtracting (7.9) from (7.8), we find that 
(b 2 - af) ^268 + b' (v 2 e - (v 2 e )) = 0. 



Hence, we may write 



¥2(0, r,t) 



b'o(r) 



o§(r) - 6g(r) 



u(9,t), 



(7.10) 



where u satisfies 
uee = v 2 e - {v 2 ). 



Using (7.6) and (7.10) in (7.7), we get 

(b'of 



(qobl /2 ip2e) T + v 6t + 



2bi /2 Lbl /2 ) TT 



v = 0. 
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Averaging this equation with respect to r, we get 
(v t + Auv e ) e + Nv = 0, 

where 



A 



(b'oY 

n 2 - h 2 
"0 °0 



dr, 



N= — 



1 I 1 



dr. 



Here, § denotes an average over a period in r. 

To make the dependence of (po on v explicit, we introduce a new time variable 
T = Wt. 



We may then rewrite (7.9) as 
V^ott + b' = 0. 



(7.11) 



Given a solution of this equation for (po(T), we have N = (vg)M, where 
1 / <?o ( 1 N 



M = 



^ bT wr, 



dT 



TT 



is a constant independent of v, and § denotes an average with respect to T 
over a period. From (7.11), we have 

\pl T + b = E 



for some constant E. Using an integration by parts, we may also write M as 



2Jh\Abl qlj 

4&§ ql 



E-b 
bo 



dT. 



Thus, the final equations for u(9,t), v(9,t) are 
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(v t + Auvg) d + M(v 2 )v = 0, (7.12) 

uee = vg- (vg). (7.13) 

It follows from these equations that, for smooth solutions, 
(vl) t + [A (uv 2 - X -u 2 + (v 2 )u e ^ + M(v 2 )v 2 ^ = 0. 

Taking the average of this equation with respect to 9, we find that {v 2 ) is 
constant in time, as stated earlier. 

In summary, the asymptotic solution of (4.3)-(4.4) is given by 



V? = Mr) + o r 5 0(T > u(9, t) + 0(e^ 2 ), 

e l/2 

^= ; 771727 - «(g,0 + Q(e), 

where the multiple-scale variables 0, r are evaluated at (7.3), s satisfies (7.4), 
(p satisfies (7.9), and u, v satisfy (7.12)— (7.13). 

If A ^ then we may rescale variables in (7.12)-(7.13) to get (7.1)-(7.2) with 
M 



As an example, let us consider the wave speeds in (4.2) arising from the one- 
dimensional director field equations, where 

b(<f) = \J f3 sin 2 ip + 7 cos 2 (p. 



If (5 < 7, then 6 has a minimum at y? = 7r/2, and (7.9) has periodic solutions for 
the mean field ip that oscillate around ir/2. Our asymptotic solution applies 
in this case. If f3 > 7, then b has a minimum at (p = 0. Although (7.9) also has 
periodic solutions in this case, there is a loss of strict hyperbolicity at tp — 0, 
where a = b, and the asymptotic solution breaks down. 
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A Algebraic details 

A.l Linearized equations 

The linearization of (2.5) at n is 

n' tt = aV(divn') - /3curl (A'n ) - 7curl (B' x n ) + A'n , (A.l) 

where 

A' — n • curl n', B' = n x curl n', 

and the Lagrange multiplier A' is chosen so that n • n' = 0. 
The Fourier mode 

n'(x, t) = n e lk - x " 1 ^, A'(x, t) = A e lk ' x "^ 

is a solution of the linearized equations (A.l) if 

£n-An = 0, n • n = 0, (A.2) 

where the linear map L : R 3 — ► R 3 is defined by 

£n = cu 2 n - a (k • n) k + (3 A (k x n ) + 7k x (B x n ) , (A. 3) 

A = n ■ (k x n) , B = n x (k x n) . 

We solve this eigenvalue problem in the next proposition. 

Proposition 5 Suppose that n , k G M 3 where n is a unit vector and k 
is non-zero, and a, (3, 7 G R are distinct positive constants. For uj G R let 
C : R 3 — > R 3 &e t/ie linear map defined in (A. 3). Then the linear system 

Cm — An = F, n • m = G (A.4) 

/10s a unique solution for {m, A} G R 3 x R /or every {F, G} G R 3 x R unless 
uj 2 = a 2 (k;n ) or u 2 = 6 2 (k;n ) ; where a 2 , b 2 are defined in (2.7)-(2.8). 

(a) Ifk is not parallel to n and uj 2 = a 2 (k;n ) ; then the general solution of 
(A./) when F = 0, G = is 

m = cR, A = -c(a-7)(k-n )(k-R), (A.5) 
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where c is an arbitrary constant and R = k — (k • no) no. Equation (A. 4) is 
solvable for {m, A} if and only if {F, G} satisfy 

R-F + (a-7)(k-n )(k-R)G = 0. (A.6) 



(b) Ifk is not parallel to n and uj 2 = 6 2 (k;n ) ; then the general solution of 
(A. 4) when F = 0, G = is 

m = cS, A = 0, (A.7) 



where c is an arbitrary constant and S = k x n . Equation (A. 4) is solvable 
for {m, A} if and only ifF satisfies 

S-F = 0. (A.8) 



(c) 7/k is parallel to n andu 2 = 7(k-n ) 2 , then the general solution of (A. 4) 
when F = 0, G = is {in. A} = {nr 1 , 0} where nr 1 is any vector orthogonal 
to n . Equation (A. 4) is solvable for {m, A} if and only ifF is parallel to n . 

Proof. First, we suppose that k is not parallel to n . Expanding 

m = TOik + m 2 R + m 3 S, 
F = F,k + F 2 n + F 3 S, 

we find, after some algebra, that (A. 4) is equivalent to 
(k • no) mi = G, 

A + (k • no) (u 2 - 7A; 2 ) m 2 = -F 2 , 
(uj 2 — ak 2 ^j mi + (uj 2 — a 2 ^j m 2 = Fi, 
(u 2 - b 2 ) m 3 = F 3 . 

The first two equations determine mi and A, and the remaining two equations 
determine m 2 , m 3 unless u 2 = a 2 or u 2 = b 2 . 

If u 2 = a 2 , then uj 2 7^ b 2 (since k is not parallel to n ) and the fourth equation 
is solvable for m 3 . The third equation is solvable for m 2 if and only if 

(k • n ) Fi - [a 2 - ak 2 ) G = 0. 

Using the equations 

a 2 - a k 2 = -(a- 7) (k • n ) 2 , R ■ F = (k • R) Fi 
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in this condition, we get (A. 6). If F = 0, G — 0, then we find that mi = 
0, m 3 = and m 2 = c, where c is an arbitrary constant. Computing the 
corresponding value of A, we get (A. 5). 

If uj 2 = b 2 , then the first three equations are solvable for A, mi, m 2 . The last 
equation is solvable for m 3 if and only if F 3 = 0, which gives (A. 8). If F = 0, 
G = then A = 0, mi = 0, m 2 = 0, and m 3 = c, which gives (A. 7). 

Finally, we suppose that k = kn is parallel to n . Then 

Cm = (uj 2 — m — (a — y)k 2 (n • m) n . 

Equation (A. 4) is therefore uniquely solvable unless to 2 = 7/c 2 , when 
Cm = —(a — -j)k 2 (no • m) n . 



In that case, (A. 4) is solvable if and only if F = Fn is parallel to n , and the 
solution is 

m = Gn + m ± , A = - \f + (a - y)k 2 G 



where m -1 is an arbitrary vector orthogonal to n . □ 

From Proposition 5, the solutions of the eigenvalue problem (A. 2) are given 
by (2.6)-(2.10). 



A. 2 Weakly nonlinear splay waves 



We look for an asymptotic solution of (2.5) of the form (3.1)-(3.2). We expand 
derivatives as 

d t ^--d e + d t , V^-d e + V, (A.9) 

e e 



where u, k are defined in (3.3). The corresponding expansions of A, A, B are 

A = -Ai + A 2 + . . . , A = Ai+ sA 2 + . . . , B = B 1 + eB 2 + . .. 

e 



where 
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A 1 = n • (k x n w ) + n • curl n , 
B x = n x (k x me) + n o x curln , 

A 2 = n • (k x n 2d ) +n r (kx me) + n • curl m + ^1 • curl n , 

B 2 = n x (k x n 2 e) + m x (k x n ie ) + n x curl m + m x curl n . 

We write these expressions as 

A i = A i + A i , Bj = Bj + Bj, (A.10) 

A = n • (k x me) , B, = n x (k x me) ■ 

We use these expansions in (2.5), Taylor expand the result with respect to e, 
and equate coefficients of powers of e. At the order e -1 , we obtain linearized 
equations for ni, Ai, which have the form 

£ni fl - Aino = 0, n ■ ni = 0, (A.ll) 
where C is the linear map defined in (A.3). 

The system (A.ll) has the splay-wave eigenvalue uj 2 = a 2 (k; no), where a is 
defined in (2.7). The corresponding solution for n 1; after two integrations with 
respect to 9, is 

ni(0,x,i) = w(#,x,t)R(x,t), 

where u is an arbitrary scalar- valued function, and R is defined in (2.9). The 
solution for Ai is 

Ai = - (a - 7) uee (k • n ) (k • R) . 

At the order e°, we obtain equations for n 2 , A 2 . Using the fact that n is a 
solution of (2.5), with A = Ao say, we may write them as 

Cmee ~ A 2 n = F 1; n • n 2 = G ± , (A. 12) 

where A 2 = A 2 — A and 

Fi = 2umet + w*nie + a {(div mo) k + V (k • me)} 

-13 {A 2e (k x no) + A 1 (k x me) + k x (A 1 n 1 ) e 

-I-Axcurlno + curl (A^q) j 
+7 {-k x (B 2 g x n ) + Bx x (k x me) - k x (Bi x n^ 
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+Bi x curln — curl (Bx x n )j + Aiiii, 

Gi = --ni • ni. 

From Proposition 5, equation (A. 12) is solvable for n 2ee and A 2 if and only if 
R • Fx + (a - 7 ) (k • n ) (k • R) G we = 0. 

After some algebra, we find that this solvability condition gives (3.4). 
A. 3 Weakly nonlinear twist waves 

We look for an asymptotic solution of (2.5) of the form (3.6)-(3.7), expanding 
derivatives as in (A. 9). The corresponding expansions of A, A, B are 

. 1 , 1. 1 , 

A = ^A 1 + A 2 + e 1/2 A 3 + ..., 
B = ^B 1 + B 2 + e 1 / 2 B 3 + ..., 
where 

Ax = n • (k x n ie ) , 
Bi = n x (k x n w ) , 

A 2 = n • (k x n 28 ) + n r (kx n ie ) + n • curl n , 
B 2 = n x (k x n 29 ) + rix x (k x n ie ) + n x curl n , 
A 3 = n • (k x n 39 ) + n r (kx n 2e ) + n 2 • (k x n w ) 

+n • curl ni + ni • curl n , 
B 3 = n x (k x n 3e ) + m x (k x n 2e ) + n 2 x (k x n 19 ) 

+n x curlnx + n x x curln . 

We write these expressions as in (A. 10). 

Using these expansions in (2.5), Taylor expanding the result, and equating 
coefficients of powers of e 1 ^ 2 , we get at the order e~ z / 2 the linearized equations 
(A. 11). From Proposition 5, this system has the twist- wave eigenvalue u 2 = 
b 2 (k; no), where b is given by (2.8). The corresponding solution for {nx, Ax} is 

ni(0, x, t) = v{9, x, i)S(x, t), Ax = 0, (A.13) 



46 



where v is an arbitrary scalar- valued function, and S is defined in (2.10). 
Equating coefficients of the order £ _1 , we obtain that 

£11200 - A 2 n = Fi, n • n 2 = G ± , (A.14) 

where £ is defined in (A. 3) and 



Fx = -p A 2g (k x n ) + k x (A 1 n 1 ) e + A 1 (kxn 



10, 



"7 



k x (b 29 x n ) + k x (Bi x ni) e - Bi x (k x n ie ) + Aiiii, 



d = --ni ■ ni. 

Using Proposition 5 to solve (A.14), we find that 



n 2 = —ru (k • n ) R + -v 2 (k x S), 
\ 2 = v 2 e {pk 2 + r(P- 1 ) (k-n ) 2 } 

where u(9,x.,t) satisfies (3.11), and 

r = . 

a — f3 



k 2 -(k-n y 



(A.15) 



We could add an arbitrary scalar multiple of the null- vector S to n 2 , but this 
would not alter our final equations, so we omit it for simplicity. 



Equating coefficients of the order e 1//2 , we obtain that 
£n 3 06» — A3n = F 2 , n • n 3 = G 2 , 



(A.16) 



where 



F 2 = 2ujTa. iet + u; t ni0 + A 2 ni + Ain 2 + a [div (n 16 ) k + V (k • n^)] 

-(3 A 3e (k x n ) + k x (A 2 n 1 ) d + k x (Ain 2 ) fl + A 2 (k x n w ) 
+A 1 (k x n 2g ) + Aicurl n + curl (Ain )] 



"7 



k x (B 3e x n ) + k x (B 2 x n 1 ) e + k x (B : x n 2 ) e 

-B 2 x (k x n ie ) - Bi x (k x n 2e ) 
— Bi x curln + curl (Bi x n )] , 



G 2 



-ni • n 2 . 
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From Proposition 5, this system is solvable if S • F2 = 0. After some algebra, 
we compute that this solvability condition gives (3.10). 



A. 4 Weakly nonlinear polarized waves 



We assume that n and k are constant and k = kn is parallel to n . We look 
for an expansion of the same form as the one used in the previous section for 
twist waves, with a plane- wave phase t) = k • x — cut. 

When k is parallel to n , the linearized dispersion relation (2.6) reduces to 
so uo 2 = 7A; 2 . From Proposition 5, the solution of the leading- 
order 0(e -3 / 2 )-equation (A. 11) is then 

ni (0,x,*) = u(0,x,*), A l = ' 

where u is an arbitrary vector such that n • u = 0. 
The 0(£ _1 )-equation (A. 14) simplifies to 



Cn 2ee ~ A 2 n = -^k 2 (u fl • u e ) n , 
n ■ n 2 = --(u- u), 

whose solution is 



n 2 = --(u • u)n , 

A 2 = i(a - j)k 2 (u • u) ee + 7/c 2 (u e ■ u e ) . 

We could add to n 2 an arbitrary vector orthogonal to n , but this would not 
alter our final equations, so we omit it for simplicity. 

Computing F 2 in the 0(e -1 / 2 )-equation (A. 16) and imposing the solvability 
condition in Proposition 5 that F 2 is parallel to n , we get after some algebra 
that 



oj _ ak 2 . . 
u et + -n • Vu 9 + — (u ■ u e ) e u 

(3k 2 



2uo 
' 2u 



{[u ■ (no x U00)] (no x u) + 2 [u • (n x u e )\ (n x u d )} 
{[(u-u)u e ] e -(u 9 - Ue )u} = 0. (A.17) 
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A computation in terms of components show that if u is orthogonal to the 
constant unit vector n , then 

[u ■ (n x u ee )] (no x u) + 2 [u • (n x u e )] (n x u e ) 

= (u • U e ) U - [(U • U) U e ] e + (ug ■ Ug) U. 

Using this result in (A. 17), we obtain (3.13). 
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